Quantum phase transitions in the Triangular-lattice Bilayer Heisenberg Model 



Rajiv R. P. Singh 

Department of Physics, University of California, Davis, California 95616 

Norbert Elstner 

Physikalisches Institut, Universitat Bonn, Nufiallee 12, D-53115 Bonn, Germany 

(February 7, 2008) 



00 

On 
ON 



3 

Hi 



> 
o 

oo 
o> 



c 

o 
o 



We study the triangular-lattice bilayer Heisenberg model with antiferromagnetic interplane coupling J± and 
nearest-neighbor intraplane coupling J = A J± , which could be ferro- or antiferromagnetic, by expansions in 
A. For negative A a phase transition is found to an ordered phase at a A c = —0.2636 ± 0.0001, which is in 
the 3D classical Heisenberg universality class. For A > 0, we find a transition at a rather large A c ~ 1.2. 
The universality class of the transition is consistent with that of Kawamura's 3D antiferromagnetic stacked 
triangular lattice. The spectral weight for the triplet excitations, at the ordering wavevector, remains finite 
at the transition, suggesting that a phase with with free spinons does not exist in this model. 

PACS: 75.10.Jm, 75.40.Cx, 75.40.Mg, 75.50.-y, 75.50.Ee 



In recent years much interest has focussed on the na- 
ture of quantum disordered phases of the Heisenberg an- 
tiferromagnets, where the combination of low dimension- 
ality, low spin and frustration, cause the ground state of 
the system to be disorderedBEl The case of one-spatial di- 
mension is relatively well studied and understood. Two 
dimensional systems have received particularly large at- 
tention, due to their relevance to high temperature su- 
perconductivity. However, despite much effort, quantum 
disordered phases are not fully understood in d = 2. 

A special situation are those types of quantum disor- 
dered phases where the ground state is to a good approxi- 
mation a product of local singlets over even-spin clusters. 
This can arise due to explicitely dimerized (or clustered) 
Hamiltonians; a situation that appears to be relevant for 
the material CaV/jOgia Another scenario is the sponta- 
neous breaking of translational symmctcry, as found in 
large- N theories and also suspected ia [Several frustrated 
models, which leads to dimerization.u □ In all these sys- 
tems, the elementary excitations are triplets with a finite 
excitation energy. 

In contrast to these, a different class of quantum disor- 
dered phases would be one where the elementary excita- 
tions are free spin-half objects or spinons. Such phases, 
for d = 2, have been predicted in systems where the 
classical ground state is non-colineaiu and their proper- 
ties have been investigated by field-theoretic methodsH. 
However, no lattice models are known where such a be- 
havior is realized. One potential candidate system for 
such a behavior is the Kagome-lattice antiferromagnet, 
where the ground-state is widely believed to be magnet- 
ically disordered. mL3 

Here we study the triangular-lattice bilayer Heisen- 
berg model. The model consists of two-layers of trian- 
gular lattices, one on top of the other, with an intralayer 
nearest neighbor Heisenberg coupling J, which could be 
ferro- or antiferromagnetic, and an antiferromagnetic in- 
ter layer nearest neighbor coupling Jj_ , between the spins 



on top of each other. Thisjaaodel maybe relevant to bi- 
layer quantum-hall systemsO and to layers of He 3 . 

The corresponding square-lattice Heisenberg-. bilayer 
has been extensively studied by many authors .113 It was 
found that with increasing A the triplet excitations at 
the ordering wavevector soften and at a critical A w 0.40 
the gap closes and there is a continuous quantum phase 
transition to a magnetically ordered phase. It should be 
mentioned that the same scenario is found for ferromag- 
netic intraplane coupling A < with a critical point at 
A c » —0.44, which, to our knowledge, has not been re- 
ported in the literature before. 

We study the triangular-bilayer model by a strong cou- 
pling expansion in the parameter A = J/Jj_. At A = 0, 
the ground state consists of products of singlets over 
pairs of spins and the elementary excitations are iso- 
lated triplets. We focus on the ordering susceptibility, the 
triplet excitation spectrum and its spectral weight as a 
function of A. For ferromagnetic intralayer couplings, the 
model exhibits a continuous quantum phase transition to 
an ordered phase at a critical coupling A c w —0.2636. 
Our results for the critical exponents associated with the 
closing of the triplet gap and the divergence of the mag- 
netic susceptibility are consistent with those of the 3D 
classical Heisenberg model. 

For the antiferromagnetic intralayer coupling, it is 
more difficult to reliably estimate the critical A c , at which 
the gap vanishes, as A c is large and the convergence of 
the series not very good at these values of A. Unbiassed 
analysis for the susceptibility and the inverse of the gap 
series shows evidence for a transition at A c w 1.2 with 
an exponent v ~ 0.53 and 7 w 1.1. These results are 
also confirmed by biassing the critical point in the se- 
ries analysis. These results are consistent with the uni- 
versality class discussed by Kawamurata for the stacked 
triangular Heisenberg antiferromagnet. 

These results also appear to rule out a phase in the 
model, at intermediate values of A, where there is still 
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a gap in the spectrum but unbound spinoas. become the 
elementary excitations. It has been shownLrB that in this 
case the order-disorder transition will lie in the univer- 
sality class of the 3-dimensional 0(4) model. One might 
also expect that in this case the triplet excitations would 
decay into the two-spinon continuum and not have a well 
defined energy momentum relation. However, we find 
that as A is increased in the model and the minimum 
of the triplet spectrum becomes more pronounced, the 
spectral weight of the triplets is reduced over much of 
the Brillouin zone but it stays finite and grows with A in 
the vicinity of the ordering wavevector. Together with 
the estimates for the critical exponents at the transition, 
this suggests that free spinons do not exist in this model. 
Our results provide further support for existence of anti- 
ferromagnetic order jp-the single-plane triangular-lattice 
antiferromagnet .E3o'E£l 

The triangular-lattice bilayer Heisenberg model is 
given by the Hamiltonian: 

H = J± Saj ■ Sba + J ^ [Sx.i • &A,j + ^B,i ■ Sb,j] 

i <i,j> 

(1) 

where A and B refer to the two layers of the triangular 
lattice and < i,j > are nearest neighbours in a given layer 
of the lattice. The triangular lattice sites are spanned by 
the two nonorthogonal primitive vectors 

ei = (1,0) and e 2 = i(-l,\/3) 

For J = 0, spins are coupled only in pairs and the 
ground state consists of product of singlets over these 
pairs. The excitations are isolated triplets localised at 
some site i. For finite values of A = J/J± an effec- 
tive hamiltonian H cS (Rj.j) describing the interaction be- 
tween these localised degenerate triplet states can be de- 
rived by a systematic expansion in A: 

H cff (R) = ^A"^„(R) 

n 

The methods for calculating H cS in powers of A =-J / J± 
are well developed and discussed in the literatureJlj The 
excitation spectrum is given by the eigenvalues of the ef- 
fective Hamiltonian W(R), where Rjj = — rj is the 
vector connecting sites i and j. It can easily be diago- 
nalised by a Fourier transform. We have calculated these 
quantities complete to 10-th order. 

In table I, the expansion coefficients for E(q) with 
q = and Qaf = 47r/3 ei are presented, correspond- 
ing to the ordering wavevectors for the ferromagnetic 
and the antiferromagnetic systems. The expansion coef- 
ficients for the.|Biagnetic susceptibilities at the same two 
wavevectors, calculated to 10-th order are given in 
the table II. 

In addition to the wavevector dependent susceptibili- 
ties and the excitation spectra, we also calculate series 
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TABLE I. Series for the triplet energies E(q) 
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TABLE II. Susceptibility series x(q) 



for the spectral weights associated with the excitations. 
The spectral weights are defined by the delta-function 
piece of the dynamical correlation function. 

S(q, w) = A(q)S(e(q) - u>) + B(q, to) 

They are calcualted via the spin-spin correlation func- 
tions, where the intermediate [States are restricted to the 
elementary triplet excitations £a These latter calculations 
are more difficult and are only done to 6th order. 

For ferromagnetic intraplane couplings, the critical 
point occurs at small values of A < 0, so even with rel- 
atively short series we can determine the critical point 
quite well and also get reasonable estimates for the criti- 
cal exponents. For the susceptibility series, the dlog Pade 
approximants lead to estimates 

A c = -0.26362 ± 0.00009, 7 = 1.407 ± 0.004. 

Here the uncertainties reflect the spread between differ- 
ent Pade estimates. Applying d-log Pade approximants 
to the inverse of the energy-gap series one obtains, 

A c = -0.2641 ±0.0005, v = 0.73 ± 0.01. 

Here, the estimates for the critical points and expo- 
nents from individual approximants are correlated and 
the more negative the critical point estimates, the larger 
is the exponent. Given the length of our series, these 
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FIG. 1. Brillouin zone of the triangular lattice 

numbers are in quite good agreement with the 3d clas- 
sical Heisenberg university class, where the best current 
estimates come from field-theory 7 = 1.3866 ± 0.0012, 
v = 0.7054 ± O.OOII.Ej That the series analysis gives 
slightly higher estimates for the exponents is common 
to man*-, mo dels and is primarily due to corrections to 
scaling E3 

We now consider the analysis for A > 0, which cor- 
responds to antiferromagnetic intraplane couplings. We 
analyzed the inverse of the energy-gap and the ordering 
susceptibility series using d-log Pade approximants and 
differential approximants. In this case, the convergence 
was much poorer as the critical point occurs at much 
larger A c . The estimates for the critical points and the 
exponents show tremendous scatter. Assuming that the 
two series have the same critical point, the most consis- 
tent estimate for A c is in the range 1.18 — 1.29. In that 
range there are four approximants for the susceptibility 
series, which give (A c , 7) values of (1.19, 1.10), (1.19, 1.08), 
(1.21, 1.06), (1.26, 1.32) and four approximants for the 
inverse gap series, which lead to (A c ,^) values of 
(1.18,0.45), (1.22,0. 51), (1.25,0. 59), (1.29,0.57). These lead 
us to conclude that A c « 1.2 and that v ss 0.53 and 
7 w 1.1. These exponents are also obtained if the approx- 
imants are biassed to have the critical point near A c = 
1.2. These results are consistent with Kawamura's uni- 
versality class for the stacked triangular-lattice Heisen- 
berg model, where he found v 1=3 0.55, and 7 « 1.1,E3 and 
not consistent with the p(4) universality class, which has 
v w 0.74 and 7 « 1. 47. M Following Chubukov, Sachdev 
and Senthil,0 such a behavior is to be expected if there 
is a direct transition from the disordered phase without 
free spinons to the 3-sublattice ordered phase. 

Another way to explore the existence of an interme- 
diate phase with free spinons is to study the spectral 
weights for the triplets and see if it vanishes as A is in- 
creased. When the spinons become the elementary exci- 
tations, the triplets can break up into a pair of spinons 
and thus will not remain sharp excitations. To analyze 
the series for the triplet spectra, we use Euler transforms 
and Pade approximants. In Fig. 1, we show the Bril- 
louin zone of the triangular lattice. In Fig. 2, the exci- 
tation spectra for A = 0.2, 0.4 and 1.0 are shown along 
selected contours. In Fig. 3, the spectral weights esti- 
mated by the [3/3] Pade are shown along the same con- 
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FIG. 2. Dispersion of the triplet excitations 
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FIG. 3. Spectral weights along selected contours 

tours. It is evident from these plots that as A is increased, 
the triplet dispersion develops a sharp minimum at the 
ordering wavevector of the triangular-lattice Heisenberg 
model. The spectral weight associated with the triplets is 
rapidly reduced over much of the Brillouin zone, however, 
in the vicinity of the ordering wavevector, the spectral- 
weights continue to increase with A, and the triplet ex- 
citations remain sharp. This provides further evidence 
for the absence of an intermediate phase in this model, 
and a direct transition from the local-singlet phase to the 
magnetically ordered phase. 

These results provide further confirmation that 
the sinek-plane triangular-lattice antiferromagnet is 
ordered.Eil However, the ratio of the ferromagnetic to 
antiferromagnetic critical points is almost an order of 
magnitude smaller for the triangular lattice than for 
the square-lattice. This, together with previous per- 
turba tiye .studies of the triangular lattice Heisenberg 
mode£j'L3 shows that the antiferromagnetic ordering in 
the triangular-lattice model is much less robust. 

In conclusion, in this paper we have studied the quan- 
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turn phase transitions in the bilayer triangular-lattice 
Heisenberg models in a strong coupling expansion. For 
ferromagnetic intralayer coupling, the transition to the 
ordered phase is found to be in the 3D classical Heisen- 
berg universality class. The antiferromagnetic intraplane 
coupling case appears to be quite different. We find ev- 
idence that there is a transition to an ordered phase at 
much larger values of A and the transition is in the uni- 
versality class of the stacked triangular lattice. This, to- 
gether with the result that the triplet spectral weight 
near the ordering wavevector continues to grow with A 
suggest that in this model, a phase with free spinons 
does not exist and there is a direct transition from the 
local singlet phase to the 3-sublattice ordered phase. This 
study lends further support to the idea that the single- 
plane spin-half triangular-lattice Heisenberg model is or- 
dered, but that this ordering is much less robust than for 
the square-lattice. 

Acknowledgements: We would like to thank Subir 
Sachdev and Olcg Starykh for valuable discussions. This 
work is supported in part by the US National Science 
Foundation under Grants No. DMR-96-16574. 



B50, 6871 (1994). 

M.P.Gelfand, Solid State Comm., 98, 11 (1996). 
H.-X. He, C. J. Hamer and J. Oitmaa, J. Phys. A 17, 1649 
(1984); M. P. Gelfand, R. R. P. Singh, and D. A. Huse, J. 
Stat. Phys. 59, 1093 (1990). 

R.R.P.Singh and M.P.Gelfand, Phys. Rev. B52, 15695 
(1992). 

J. C. Le Guillou and J. Zinn- Justin, Phys. Rev. Lett. 39, 
95 (1977). 

A. J. Liu and M. E. Fisher, Physica A 156, 35 (1989). 
P. Azaria, B. Delamott and T. Jolicoeur, Phys. Rev. Lett. 
64, 3175 (1990). 



1 S. Chakravarty, B. I. Halperin and D. R. Nelson, Phys. 
Rev. B39, 2344 (1989). 

2 A. V. Chubukov, S. Sachdev and J. Ye, Phys. Rev. B 49, 
11919 (1994). 

3 K. Ueda et al., Phys. Rev. Lett. 76, 1932 (1996); M. P. 
Gelfand et al., Phys. Rev. Lett. 77, 2794 (1996); W. E. 
Pickett, Phys. Rev. Lett. 79, 1746 (1997). 

4 N. Read and S. Sachdev, Phys. Rev. Lett. 62, 1694 (1989). 

5 M.P.Gelfand, R.R.P.Singh and DA. Huse, Phys. Rev. B40, 
10801 (1989). 

6 H.J.Schulz and T.A.L.Ziman, Europhysics Letters 18, 355 
(1992). 

7 A.V. Chubukov, S. Sachdev and T. Senthil, Nuc. Phys. 
B426, 601 (1994). 

8 A.V. Chubukov and OA. Starykh, Phys. Rev. B53, 14729 
(1996). 

9 P. Lecheminant et al. Phys. Rev. B56, 2521 (1997). 

10 R.R.P.Singh and DA. Huse, Phys. Rev. Lett. 68, 1766 
(1992). 

11 S. DasSarma, S. Sachdev, and L. Zheng, Phys. Rev. Lett. 
79, 917 (1997). 

12 K. Hida, J. Phys. Soc. Jpn. 61, 1013 (1992); A.W. Sandvik 
and D.J. Scalapino, Phys. Rev. Lett. 72, 2777 (1993); M.P 
Gelfand Phys. Rev. B53, 11309 (1996); W.H. Zheng, Phys. 
Rev. B55, 12267 (1997). 

13 H. Kawamura, J. Phys. Soc. Jpn. 61, 1299 (1992); ibid. 59, 
2305 (1990). 

14 B. Bernu, P. Lecheminant, C. Lhuillier and L. Pierre, Phys. 
Rev. B 50, 10048 (1994). 

15 N. Elstner, R.R.P.Singh and A.P.Young, Phys. Rev. Lett., 
71, 1629 (1993); N. Elstner and A.P.Young, Phys. Rev. 



4 



